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ABSTRACT 

The Sunyaev Zel'dovich (SZ) signal is highly non-Gaussian, so the SZ power spectrum (along with 
the mean y parameter) does not provide a complete description of the SZ effect. Therefore, SZ-based 
constraints on cosmological parameters and on cluster gastrophysics which assume Gaussianity will 
be biased. We derive an analytic expression for the n-point joint PDF of the SZ power spectrum. Our 
derivation, which is based on the halo model, has several advantages: it is expressed in an integral 
form which allows quick computation; it is applicable to any given survey and any given angular scale; 
it is straightforward to incorporate many of the complexities which arise when modeling the SZ signal. 
To illustrate, we use our expression to estimate p(Ce), the one-point PDF of the SZ power spectrum. 
For small sky coverage (applicable to BIMA/CBI and the Sunyaev Zel'dovich Array experiments), our 
analysis shows that p(C£) on the several arc-minute scale is expected to be strongly skewed, peaking 
at a value well below the mean and with a long tail which extends to tail high Ct values. In the 
limit of large sky coverage (applicable to the South Pole Telescope and Planck), p(Ct) approaches a 
Gaussian form. However, even in this limit, the variance of the power spectrum is very different from 
the naive Gaussian-based estimate. This is because different I models are strongly correlated, making 
the cosmic variance of the SZ band-power much larger than the naive estimate. Our analysis should 
also be useful for modeling the PDF of the power spectrum induced by gravitational lensing at large 
t. 

Subject headings: Cosmology: theory-large scale structure-cosmic microwave background- methods: 
statistical 



1. INTRODUCTION 

The Sunyaev Zel'do vic h (SZ) effect 
(IZeldovich fc Sunvaevl Il969t ISunvaev fc Zeldovichl 
Gaza is a powerf ul probe of baryons, dark matter 
and dark energy (|Birkinshawl Il999t ICarlstrom et al.l 
2002). To extract the full statistical power of this 
probe, SZ statistics must be understood to an accuracy 
which matches that which the observations will reach. 
The statistic which has been most actively studied 
is the ensemble average power spectrum Ci of the 
thermal SZ signal. For the simp lest case of ad i abatic 
gas evolution both analytic (ICole fc Kaiser) 1 19881 : 
Makino fc Sutol Il993t l Atrio-Barandela fc Mucketl Il999i 



Komatsu fc Kitavamal 1999; 



2000t iMaiumdarl [200lT 



Molnar fc Birkinshawl 
IZhang fc PenT 



2001; 



Komatsu fc Seliakl l2002f) and numerical (Ida Silva et 
20001 iR cfrcgicr ct al. 2000; Seli ak Burwell fc Penll200lL 
Springel. White fc Hern auist 200ll: IZhang. Pen fc Wang 
2002T) approaches are beginning to give consistent res ult s 



al. 

31; 



( Refr egier fc TevssieHl2002l: IKomatsu fc Seliakll2002D . 

However, Ct alone (or in combination with the 
mean y parameter) does not provide a complete de- 
scription of the SZ effect. This i s becau s e the 



SZ 



is 



non-Gaussian (ICoorav 
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: lYoshida. Sheth & Diaferio 


20011 IZhang. Pen & Wang 
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2002). In contrast to the primary CMB and weak grav- 
itational lensing sig nals, the SZ signal is non-Gaussian 
even on large scales (I Zhang fc Pen|[200ll ). This is because 
the effect is dominated by rare massive clusters. There- 
fore, interpretation of a SZ Ct measurement requires 
not only a_ robust prediction of the expected ensemble- 
averaged Ce, but also the probability distribution func- 
tion (PDF) of Ct around this mean value. 

Because the SZ effect is dominated by the most mas- 
sive clusters, and the abundance of such clusters is 
very sensitive to the am p litude of density fluctu ations 
([Press fc Schechterl Il974t ISheth fc Tormenl 1 1999( 1. the 
SZ power spectrum depends strongly on a$: Ci oc 
<7g 7 for the range of currently interesting cosmological 
models dSeliak. ButwcU TTctI2001I : [Zha7g fc Pedl200ll : 
IKomatsu fc Seliakl l2002f ). In principle, this makes the 
SZ power spectrum a powerful probe of a ^. The CMB 
powe r spectra excess repo rted by CBI dMason et al. | 
' 20031: iReadhead et all 12004( 1 and BIMA ([Dawson et all 
2002T ) are consistent with having a contribution from 
the SZ effect provided cr 8 ~ 1.0 dDawson et all 120021 : 
iBond et alJ l2005t IReadhead et~all |2004j). ThisSZ ex- 
planation is further supported b y the frequency depen- 
dence of the CMB power excess ([Kuo et al.ll2006l l. How- 
ever, this value is uncomfo rtably high compar ed to the 
WMAP3 value of a s « 0.8 (iSnergel et al.ll20T)l . 

The SZ-based estimates of as a ssume adiabat i c evo- 
lution of the gas. Feedback (|da Silva et alJ 120011 : 



IWhite. Hernquist fe Spingell 120021: iLin et al.l I2004D and 



cooling ( da Silva et al.1 l200ll iZhang fc WJ T2003) effects 
can decrease the SZ power spectrum by a factor of 2; 
if these are included, then the required value of as in- 
crease s further. O ther mechanisms such as magnetic 
fields ljZhangll2004fl can further reduce the amplitude of 
Cj. Other contributions to the SZ sign al have also bee n 
studied— the SZ effect of the first stars (|Oh et al.ll2003l ). 
and contributions from unvirialize d intergalactic medium 
(|Atrio-Barandela fc Mucketl |2006|) . But if the WMAP 3 
estimates are accurate, r ~ 0.1 and as — 0.8, then these 
additional contributions are likely to be subdominant. 
Thus to explain the CBI/BIMA SZ measurements it ap- 
pears that a 8 > 1.1-1.2: evidently, a severe discrepancy 
with the CMB-based estimate exists . 

There are at least two ways in which the discrepancy 
can be decreased. One is to assume that the uncertainty 
on the estimate is artificially small. This would happen, 
for instance, if the statistical errors on the SZ Ct and 
Cgi are highly correlated, even for widely separated £, £ 
pairs. If so, then the cosmic variance of the measured 
band-power SZ Cg would be much larger than the naive 
Gaussian-based estimates. 

Alternatively, the SZ-based estimates of as arise from 
requiring the measured Cg to match the ensemble aver- 
aged power spectrum predicted from theory. However, if 
the distribution of Cg is skewed, then this requirement 
is unreasonable. Since we have good reason to expect 
the signal to be non-Gaussian, it may be that the dis- 
crepancy can be resolved if one accounts for this non- 
Gaussianity. For small sky coverage in particular, p{Cg) 
can be highly skewed, with a non-negligible tail toward 
high C;. Thus, for a survey with small sky coverage, the 
probability of obtainin g Cg several t i mes la rger the mean 
Cg is not negligible. iDawson et alj (|2006| ) simulated an 
SZ map for one choice of cosmological model and then 
scaled the simulated SZ decrement (assuming the signal 

7 /2 7 

scales as a 8 , so that Cg scales as Og) to estimate the 
effect of this non-Gaussianity as as varies. They used 
this procedure to argue that as ~ 0.7 may well produce 
a long enough tail of Cg values to explain the CBI/BIMA 
measurement. While this is reassuring, it is not obvious 
that this simple rescaling is indeed an accurate descrip- 
tion of how the non-Gaussian distribution of Cg depends 
on as- 

The discussion above shows why understanding the 
full PDF of Cg is crucial. Current theoretical under- 
standing of p(Cg) is limited. Even for a single cosmo- 
logical model and a given sky coverage (e.g. 1 deg 2 or 
larger), simulations lack sufficient realizations to mea- 
sure the full PDF reliably. The scaling method of 
IDawson et al.l (|2006T ) only provides a crude estimate, 
since the non-Gaussianity (and the PDF) is certainly 
a function of as- Analytical estimates of the low or- 
der moments of p(Cg), such as the variance and covari- 
ance of Cg have been made (ICooravll2001t IZhang fc Pen! 
I2001 IKomatsu fc Seliakl [2002). These calculations pro- 
vide quantitative estimates of the errors in Cg, and repre- 
sent important first steps towards understanding the SZ 

Contaminations of unresolved and unremoved point sources 
remain a possible solut i on to the power e x cess problem 
IWhite fc Maiumdarl l200l IToffolatti eTaLl [20051; IDouspis et~aTl 
2006), although current modeling is quite uncertain. 



non-Gaussianity. However, to specify p(Cg), higher or- 
der moments are also required. The computational time 
required to calculate higher order moments rapidly be- 
comes prohibitive. To estimate p(Cg) for any given cos- 
mology and given survey strategy in a reasonable amount 
of time, going beyond this moment by moment calcula- 
tion is essential. 

This paper presents a fast, complete, analytical 
method for calcul ating the full n-point P DF. It is based 
on the halo model (|Coorav fc She th 2002) and allows one 
to easily estimate the effect on the SZ signal of different 
treatments of the relevant gas physics. <J2] outlines the 
halo model calculation of the lowest order statistics of 
the SZ power spectrum: the mean, the variance and the 
covariance. An analytical expression for the n-point joint 
PDF of the SZ power spectrum is derived in We show 
explicitly that it correctly reproduces the usual expres- 
sions for the mean, variance and covariance. We then 
show that the analysis is particularly simple in the limits 
of very large and very small sky coverage. A numerical 
calculation of the one-point PDF, and a comparison with 
the large sky coverage limit is presented in £13.61 £0] sum- 
marizes our results, highlights several key simplifications 
in our approach, and discusses possible extensions. 

Where necessary in this paper, we adopt a flat ACDM 
cosmology with Q m = 0.3, Oa = 0.7, a 8 — 0.9, and 
H = lOO/i km s" 1 Mpc -1 with h = 0.7. To illustrate 
our results we use an initial power spectrum with index 
n = 1, the BBKS fittin g formula for the transfer func- 
tion, the NFW profile (jNavarro. Frenk fc White! Il996h 
for the dark matter profi le, the Press-Schechter formula 
(|Press fc Schechteill 19741) for t he halo mass functio n and 
the associated halo bias factor (|Mo fc Whitdll99l when 
modeling halo clustering. When our formalism is used 
to interpret simulations or observations, more accurate 
m odels of halo abundance and clustering, such as those 
of lSheth fc Tormenl (|1999t ). should be used. 

2. THE MEAN, VARIANCE AND COVARIANCE OF THE SZ 
POWER SPECTRUM 

The dominant contributi on to the SZ power spectrum 
is from virialized regions (ISeliak. Burwell fc Pen] 1 2001 : 
White. Hernquist fc Spinge]||2002tlCoorav fc Shethl l2002 



Hernandez-Monteagudo et al.l " l2006h . This makes the 
halo model particularly well suited for estimating Cg and 
its distribution. On the angular scales of interest in what 
follows, the two-halo term is always much smaller than 
the one-h alo term: the contribution i s less than 1% at 
£ > 1000 (jKomatsu fc Kitavamalll999t ). so it can be ne- 
glected. Thus, the measured power spectrum C = C(£)\s 
at multipole £, is just a sum over the contributions from 
each cluster in the survey. I.e., for a survey with sky 
coverage / sky , 

C\s = f£ / dz^ J dM-^S(M,z;X)p(X\M,z)dX . 

(1) 

Here, V is the survey volume. The explicit factor / sk y 
in the above equation cancels the implicit dependence of 
V on / s ky (V cx /sky) so the expectation value of C\s 
does not depend on / s ky S(M,z;X) is the SZ power 
spectrum of a single cluster at the corresponding £. In 
this and the following expressions, we assume that no 
confusion will arise from the fact that we do not write 



3 



explicitly that both S and C depend on i. The quantity 
S describes the structure of the cluster; it is mainly de- 
termined by the cluster mass M and redshift z, but it 
also depends on other parameters A, which may include 
the cluster concentration parameter c, cluster shape, etc. 
These parameters may or may not correlate tightly with 
M, and, in general, they almost certainly are not fixed 
by M alone. We use p(X\M, z) to denote the distribution 
of these parameters at fixed M and z; it is normalized so 
that J p(X\M, z) dA = 1. For simplicity, we neglect these 
extra parameters in the numerical calculations which fol- 
low, although we do include them in our analytical ex- 
pressions. We will discuss their possible effects later. 

If dn/dM is the ensemble averaged halo mass function, 
then the comoving number density of clusters per mass 
interval that are actually in the survey volume is 



dn/dM = dn/dM l + S L (M)+6 



(2) 



Here, 5l(M) denotes the mean fluctuation in cluster 
number density which arises from the fact that the den- 
sity of dark matter in the survey volume may not be 
exactly the same as the universal average, and Sp rep- 
resents the fluctuation which arises from the fact that 
clusters represent a discrete point process realization of 
the smooth density field dn/dM [1 + 6l{M)]. It is stan- 
dard to assume that Sp is drawn fr om a Poisson distribu- 
tion with mean dn/dM \ 1+6l(M)] (|Casas-Miranda et al.l 
[200llHu fc Cohijl2006tlHolderll200^ ~ 

In what follows, we will often consider bins in redshift, 
mass and A. We use the Latin letter j to denote jth 
redshift bin and the Greek letter a to denote the ath bin 
of M and A. The expectation value of C\s is 



C = 



dz — — 

dz 



f dn 

dM—S{M,z;X)p{X\M,z)dX, (3) 



where V? is the comoving volume of the survey. The 
variance of C\s is the sum of two terms: 

r rjvf f dfi 
a c = f^ I dz ^ dM—S 2 (M,z;X)p(X\M)dX 



C J sky 



dM—b(M,z 3 )S{M,z f ,X) 



x p(X\M, Zj )dX AVf 



(4) 



The first term on the right hand side is often called the 
Poisson term. If 



S 2 pdndV f dX 



SpdndV 1 dX 



(5) 



denotes an effective number of clusters, then the frac- 
tional fluctuation in Ct induced by the Poisson term is 

The second term on the right hand side of equation (01 
accounts for cluster clustering. The sum is over all red- 
shift bins j. We choose bin sizes Az ~ 0.1 so that 6s, the 
mass overdensity fluctuation smoothed over the survey 
volume f a kyAVf in any given redshift bin, is (approxi- 
mately) uncorrelated with 6s in other bins. The quantity 



6 2 s .) = f P j (k)W 2 (k)d 3 k/(2n) s 
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Fig. 1. — The variance of the SZ power spectrum. The lower 
panel shows the weighted cluster number Np. (Np / s ky) — 1 ' 2 is 
the fractional statistical error caused by the Poisson distribution 
of clusters. The upper panel shows the ratio between the statistical 
error caused by the large scale clustering and the Poisson distri- 
bution. We show below that when Np f s ^ y 3> 10, then p(Ci) will 
approach a Gaussian form, in accordance with the central limit 
theorem. 



is the rms of the matter over-density when smoothed over 
the survey window in the jth redshift bin; Pj and Wj 
represent the 3D matter density power spectrum, and the 
Fourier transform of the survey window function, both 
in the jth redshift bin. The factor b represents the fact 
that the cluster distribution is biased relative to the dark 
matter: 6 L (M) = b(M)6 s . 

We model S(M,z) usi ng the universal gas profile of 
(|Komatsu fc Se liak 200lJ). This model makes three rea- 
sonable approximations/assumptions: (1) Intra-cluster 
gas is in hydrostatic equilibrium; (2) the gas pressure 
P g oc p^ and the polytropic index 7 is a constant; (3) 
the ratio of the gas density p g and the dark matter den- 
sity approaches the cosmological value Qb/QnM in the 
outer regions of clusters. Given a dark matter profile, 
this model predicts S uniquely. When combined with 
a model for the m ass function dn/dM, this model com- 
pletely specifies C l|Komatsu fc Seliak 2002T). 

Fig. Q] shows the Poisson and halo clustering contri- 
butions to the SZ power spectrum variance for a range 
of choices of the fraction of sky covered / s ky The top 
panel shows that the halo clustering contribution is sub- 
dominant, but non-negligible, at the relevant scales. Al- 
though the Poisson term has a simple dependence on 
the sky coverage (oc f~^ y ), the halo clustering term is 
more complicated. It depends implicitly on / s k y because 
a 2 depends on the shape of the survey volume. Since 
the effective power index of CDM power spectrum n e g 
varies with scale, a 2 has non-trivial dependence on the 
smoothing scale. Thus the ratio between the halo clus- 
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Fig. 2. — The cross correlation coefficient r between different 
I modes. The dependence of r on sky coverage is negligible; we 
have assumed coverage of 200 deg 2 . The value of £ for each line 
is indicated by the lowest I value for which it is plotted. 



tering term and the Poisson term has a complicated de- 
pendence On /sky 

The bottom panel shows that Np , the effective number 
of clusters, increases with increasing £. This reflects the 
fact that the small scal e SZ power spectrum is do minated 
by smaller halos (e.g. iKomatsu fc Se liak (|2002f) ). Since 

— 1/2 

the Poisson contribution to the power scales as -/V p , 
it decreases as I increases. Because the bias factor also 
decreases with decreasing halo mass, the halo clustering 
contribution also decreases on small scales. However, the 
decrease as not as dramatic as for the Poisson term; this 
is why the ratio of this to the Poisson term increases with 
I. 

The universal gas model is spherically symmetric, so 
different £ with the same £ = \£\ are identical. 2 Hence, 
the combination of {2£ + 1) modes all having the same \£\ 
has the same statistical error ac and the fractional error 
— [^Vp/sky] 1 ^ 2 - For a Gaussian random field, the frac- 
tional statistical error of the power spectrum for each I is 
[(2^ + l)/ sky /2]- 1 /2. since N P > £ (Fig. [J), the statisti- 
cal error of the SZ effect of each I is smaller than it would 
be in a Gaussian random field. The halo model shows 
that this is because clusters are nearly self similar, so the 
statistical error in the SZ power spectrum is mainly de- 
termined by the fluctuations in cluster number. Because 
there are many clusters across the sky, these fluctuations 
are small, and the resulting statistical error is smaller 
than the corresponding Gaussian estimate. 

However, this does not mean that the statistical preci- 
sion of the SZ measurement exceeds the Gaussian limit. 
Although clusters are nearly self similar, the contribu- 

2 Accounting for cluster ellipticities breaks this degeneracy as 
we discuss in [|3] 



tion of each cluster spans a large range in I. As a re- 
sult, the observed signal in the different I modes can be 
highly correlated. If A£ c denotes the correlation length 
in £ space, then it is the combined error on A£ c modes 
which is the same as that of one single mode. This can 
be much larger than the corresponding Gaussian error 
[(2£+l)A£ c W2]- 1 / 2 . 

To estimate A£ c we compute the covariance between 
CV and C»< : 



Cov„ 



/sky / dz 



(ill I 

dM — SS p(X\M)d\ 



dV f 

dz 



dM^-bSp{\) d\ AVf 
dM Hy ' 3 



dM-^bS'p(X) dXAV/ 



(7) 



where S = S e ' (M, z; A). Again, the Poisson term dom- 
inates this covariance. The cross correlation coefficient 
between the statistical errors of Ce and C'g is 



„2 _ 



o- 2 c o- 2 c , 



(8) 



Since the Poisson term dominates in both the variance 
and covariance term, the dependence of r on / s k y is neg- 
ligible. 

Fig. [5]shows r for 200 deg 2 sky coverage. Note that the 
correlation length A£ c > 1000. In simulations, the finite 
size of the SZ map map , makes the natural bin size A£ = 
2n/6 map ~ 300 for typical map ~ 1°. Thus, £ modes over 
this A£ can be treated as completely correlated. The 
fractional statistical error on C measured in simulations, 
with respect to the Gaussian random field, is expected 
to be ~ [£A£/N P ] 1 / 2 . This is indeed what is seen. For 
example, the anal ysis above suggests that the quantity 
shown in Fig. 6 of Zh ang. Pen fc Wangl (|2002j ) is 04 ~ 
{21 + l)A£[l + LSS/Poisson]/Afp - 2 (if the contribution 
of halo clustering has been taken into account). The 
numbers from Fig. Q] suggest that 04 ~ 50 at £ = 10 3 
and 04 ~ 2 at £ = 10 4 , co nsistent with the simulations 
iZhang. Pen fc Wand (pOOl) . 

3. THE PDF OF THE SZ POWER SPECTRUM 

The quantities C, ac and Cov u > are the lowest order 
moments of the power spectrum PDF p{Ce). However, 
because the SZ signal is non-Gaussian, these low order 
moments do not uniquely determine p{Cg). Since p{Ct) is 
required to make unbiased error analyses and parameter 
constraints, this section provides an analytical integral 
formula for p(Cg). This integral form, which is based on 
the halo model, allows fast numerical calculation and has 
the flexibility to incorporate much of the relevant physics. 
For clarity, § 13.11 considers the case of a single redshift 
bin; the contributions from different bins are summed in 



3.1. Single redshift bin 

Given a survey and a redshift bin, the survey volume is 
fixed. We sort the clusters in this volume into bins of M 
and, in principle, bins in the other parameters A/. Then 
the overall SZ power spectrum is 



C^/skyE^^ 



(9) 
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where N a is the number of clusters in the a-th bin and 
S a is the SZ power spectrum of each cluster in the a-th 
mass (and A) bin. The PDF of C in the jth redshift bin 
is 



Pi 



(C)= E E P(Ni,---,N a ) 



JVi=0 N a =0 



(10) 



(11) 



where P(N\, . . . , N a ) is the joint distribution of the N a . 

This distribution is computed in two steps. First, if 
the matter overdensity 8 in the survey volume is not ex- 
actly ze ro, then the mean number of expected clust ers is 
altered dMo fc Whitdfl99^ ISheth fc Tormedl2002T) . We 
use 

N a , L = N a (l + 8 a ) (12) 

to denote this altered mean. We then convert the smooth 
field N 0: l to a discrete point process by assuming a Pois- 
son sampling model. Namely, we assume that given 8 a , 



Pp(N a \d a ) 



NJ 



■exp[— JV a i]; 



(13) 



note that this Poisson process in the ath bin is indepen- 
dent of that in the other bins. Hence, 



P(N! 



,N a ) = 



Y[P P (N a \S a ) 

a 

x Pl(Si, - ■ ■ ,S a )Y[d5 a 



(14) 



To proceed, we need a model for Pl{8i, ■ ■ ■ , 8 a ), the joint 
PDFof(5 Q . 

When the volume is small, stochastic discreteness ef- 
fects dominate, and Pl(8i, • • • , 8 a ) — * J\ a Ph{S a ). When 
this limit applies, the Poisson fluctuation is much larger 
than the matter fluctuation and the exact form of Pl 
is not necessary. However , Pl(8i, ■ • ■ ,8 a ) in the general 
case is more complicated (|Sheth fe Lemsonlll999h . 

When 5 a -c 1, then the linear bias model 

8a = b a 8 (15) 

works well (IMo fc White! fl99& ISheth fc Tormenl 120021) . 
and 

(16) 



Pl ^PsJ{8v{8 a -b a 8) 



Typically, this requires the volume to be larger than 
about (100/i _1 Mpc) 3 . In practice we will apply the linear 
bias model to all relevant cases. 

To see what linear bias implies, note that the charac- 
teristic function is 



Pj ; 



dC Pj (C)exp(ifC) 
I PL(8i,---,S a )Y[dS a 

a 

E^x(e i/S — i) 

a 

d8P,j{8) cxpj / |exp(i/S(A) - 1 dN L 



(17) 



x exp 



where / = f / f s ky The second equality uses the Poisson 
model, and the hnal expression uses the linear bias model 
in the limit that the mass (and A) bin size goes to zero. 
I.e., dNi, = dN(l + b8) is the expected number of clusters 
in the survey volume and redshift and mass bin (with- 
out Poisson fluctuation), and dN / f s ky is the expected 
number of clusters averaged over the whole sky. Our no- 
tation makes explicit that the distribution Pj(8) of the 
mass density fluctuations 8 may depend on redshift. 

On large scales, the distribution of 8 should be ap- 
proximately Gaussian: Pj(8) <x exp(—8 2 /2a^). Hence, 
for each redshift bin, 



Pj=exp /sky-D, 



f 2 

J sky 



h 3 3 



where 



and 



D , 



sky 



dN (e 



ifS 



E 3 = /sky I dNb (e ifs 



1) 



(18) 



(19) 



(20) 



Notice that neither Dj nor Ej have any dependence on 
/sky, because the implicit dependence of N on / s k y (N oc 
/sky) cancels the factor in the definitions of Dj and 
Ej . We will soon show that Dj describes the discreteness 
of the cluster distribution and Ej describes the clustering 
effect. 

Strictly speaking, equation (|18[) follows from integrat- 
ing the Gaussian from — oo to +oo. Since b8 < —1 is 
problematic (one may not have a negative number for 
the mean halo count in a cell), one might have thought 
the range of integration should be restricted to 8 > —1/6. 
In practice, when the survey area is large, then Oj <C 1, 
so 8 <C 1 and bd > — 1 almost surely, and ignoring the 
restriction on 8 is reasonable. This approximation sim- 
plifies the calculation in the next subsection significantly. 

Before moving on, suppose that S was the same con- 
stant for all clusters. In this case, 



fskyDj 
j sky Ej - 



5-1 dN = [S -UK 



5-1 



dNb: 



Voir 



5 - 1 ALff 6, 



ff Ocff, 



where we have defined 5 = e^ s . Hence, 
p — » exp 



(S-l)N eS +^—^-m.bi a a 2 



eff u eS°j 



(21) 



(22) 



this expression is essentially the generating function of 
unweighted counts in cel ls. Notice that it agre es with 
the ex pressions derived bv lHu fc Cotml (f2006ft and|Holder] 
(|2006h (although Hu & Cohn set the lower limit of the 
integral over 8 to — 1/6). In this respect, our analysis 
generalizes their work to the case in which clusters con- 
tribute different weights. 

In the counts in cells case, suppose that each cluster 
is weighted by the number of galaxies it contains, and 
that W(S) is the generating function of galaxy counts per 
cluster. Then the generating function of galaxy counts in 
cells, say g(S), is given by g(S) = G(W(S)), where G{S) 
is the generating function of unweighted cluster counts 
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(e.g. our Eq. [22]). We have used W(S) to emphasize that 
one can think of the galaxy counts in cells distribution 
as arising from a process where each cluster is weighted 
differently (i.e., by the number of galaxies in it). 

In our SZ calculation, each cluster has a different 
weight, so we can think of W(S) as having a different 
value for each cluster. So our Eq. [T7]is indeed equivalent 
to the counts in cells analogy. This analogy also shows 
how to generalize our formalism to include scatter in halo 
concentrations — > scatter in SZ signal, — > scatter in SZ 
Cg even at fixed halo mass. 

One might have thought that equation (fT8| fails in the 
limit of small sky coverage. Appendix [A] shows that, in 
fact, it agrees with the Poisson (small sky coverage) limit 
to order a 2 . This agreement gives us the confidence to 
apply equation (|I8p to any relevant survey area. 

3.2. Projection along the line of sight 

The SZ power spectrum is obtained by summing over 
all redshift bins. Since the contributions from differ- 
ent redshift bins are correlated, the relation between the 
PDF of the total SZ power spectrum and the PDF of 
the SZ power spectrum from each redshift redshift bin is 
complicated. However, for a sufficiently wide bin in z, 
e.g. Az = 0.1, the fluctuations S z of different redshift 
bins are approximately uncorrelated. In this case, the 
Fourier transform of the SZ power spectrum PDF is just 
the product of equation (fT8|) for each redshift bin, 



P 



:exp 



[<?(/)] ■ 



7712^2 f2 " 



(23) 



Here, j denotes the j-th redshift bin. 

The Inverse Fourier transform of equation 
the PDF 



P(C t 



Pexp(-ifC i )df/(2n). 



yields 



(24) 



Equation (f2"3")l has all desired quantities: (1) P*(f) = 
P(-f), so the PDF is real; (2) P(f = 0) = 1, so 
/ P(C)dC = 1, meaning the PDF is correctly normal- 
ized. It is straightforward to verify that dP/df\f=Q = iC 
and d 2 P/df 2 = — ((7 2 +0-^) reduce to the expressions for 
C and Op given in the previous section (equations [3] and 
^respectively). This shows explicitly that our expression 
for P(C) has the correct first and second moments. 

3.3. 2-point joint pdf of Ci and C e > 

The joint pdf of P 2 (Ci,Cfi) can be derived similarly. 
The 2D Fourier transform is 



P 2 (fJ)= I I P 2 (C,C)exp[ifC + if C ) dfdf . 

(25) 



For a single redshift, 

P 2 = fPh(Si,---,8 a ) H d5 a 



(26) 



x exp 



£A^(e i/X+i/ ' S «--l 



Again, assuming the linear bias model and Gaussian 
Pj(5) yields 



lnP 2 = E f sky D 2 j 

j 

where 



f2 p2 J2 
hky-^2] °~j _ 



G 2 (/,/), (27) 



D 2 = L 



and 



sky 



E2 - / Sk y 



dN 



dNb 



JfS+if s _ 1 



,ifS+if s _ 1 



(28) 



(29) 



It is straightforward to verify that this expression 
is normalized to unity, and that the expressions for 
C, C ,crc, o~c' ano - Gov U ' which it implies are all in agree- 
ment with our previous expressions (e.g. equation EJ). 

If G 2 (f, f ) can be written as the sum of two terms, one 
a function of / and the other of / , then this would imply 
that C and C are independent. The expression above 
shows that, in general, this condition is not satisfied. 
For example, there are cross terms in both D and E 2 
which imply correlations between C and C . In the limit 
that t and I are sufficiently different that S ^> S (or 

S <C S ) , then D 2 roughly meets this condition. Since the 
correlation induced by the E 2 term is weak, the resulting 

P 2 {C,C')~P{C)P{C'). 

3.4. n-point joint pdf of Ce and C t > 

The analysis above for the Fourier Transform of the 
2-point pdf is easily generalized to that for the n-point 
pdf: 



lnP n (/i, • • • , f n )=ln 



x 

E 



fskyD 



f2 p2 2 
Jsky-^nj °j 



df n 



"J 



where 



Dn = t 



sky 



and 



En = /, 



1 

sky 



dN 



dNb 



,i/iSi+-+i/„S„ 



- 1 



3 i/iSi + -+i/„S„ 



1 



(30) 



(31) 



(32) 



Equation (|30|) for the 71-point PDF and its special cases, 
equations (|23[) and (|2T|) for the 1- and 2-point PDFs, are 
the key results of this paper. 

3.5. The large sky coverage limit 

To better understand the physical meaning of the 
terms in equations ([23]) . (|2T[> and (|30)) . it is useful to 
study the limit of large sky coverage. 

Consider the one point pdf P(C) in this limit. Ef- 
fectively, only those modes with Re[G(/)] > —10 con- 
tribute to P. Over this range of /, Re[G] decreases 
with /. Thus, those modes with Re[G] > —10 have 
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Fig. 3. — The requirement for having a Gaussian PDF. Solid 
lines show Re(J3 ■ Dj) for several multipoles. Straight dotted lines 

are <x / 2 . The deviation of the dotted lines from the corresponding 
/ modes implies that the pdf is non-Gaussian. The / modes below 
the straight dashed lines are the relevant modes for corresponding 
sky coverage (SPT: ~ 4000 deg 2 ; ACT/APEX: ~ 200 deg 2 ; SZA: 
~ 10 deg 2 and CBI: ~ 1 deg 2 ). 



/ < f c , where Re[G(/ c )] = -10. Since |Re(G)| in- 
creases with /sky, f c decreases with / s k y - Thus, for suffi- 
ciently large / s ky, fc and so all relevant / are small, and 
exp(i/G) - 1 ~ i/C - (/G) 2 /2. In this case, 



G(/)~i/G- 
=i/G- 



f 2 c 2 

Re(Dj) ~ -i — J-, 
y 3J 2 TV/ 



0, 

1 



and 



Nvjfsky + lj<7j . 



Pol 



Here, 



Ci=f~^ J S(X)dN J , N Pj 

, _ jS(X)b(X)dNj 
J " / 5(A) dNj 



x [fS(X)dN j ] 
/sky J S(X) 2 dNj 



(33) 



(34) 



and 



(35) 



are all defined in jthe j-th redshift bin. G is defined by 
equation ([3]) and C — J^Cj. None of these quantities de- 
pends on / s ky. Notice that the coefficient of / 2 /2 in equa- 
tion (f3"4"| is just the variance of G (equation 0]) . Thus, 
in the large sky limit, there arc sufficiently many clus- 
ters that the Poisson distribution becomes Gaussian, in 
accordance with the central limit theorem. Thus P, the 
Fourier transform of exp[G], will also be Gaussian. 



When the above conditions are satisfied, the Fourier 
transforms of equations (|23|) and (|34|) yield 
1 ( (G - G) 2 ' 



P(C) 



exp 



2a% 



(36) 



where G and ac are given by equations © and ([!]), re- 
spectively. This shows that, in the limit of large sky 
coverage, P(C) becomes Gaussian. 

However, the conditions (|3"3"]) are not always satisfied. 
The condition Re(Dj) ~ (/// s k y ) 2 G 2 /A^- is especially 
hard to meet. Since the term Re(J^ • Dj) dominates that 
involving J2j Ej (consistent with Fig. [1]), requiring that 

^Re(Z) J )oc/ 2 when 



Uy^MDj) > -10, 



(37) 

is a suitable condition for Gaussianity. We call this the 
Gaussian condition. Fig. [3] shows that, for an all sky 
survey, all I > 500 satisfy the Gaussian condition given 
above. For SPT, the Gaussian condition is satisfied at 
£ > 10 3 . However, for small sky coverage survey such as 
SZA, the Gaussian condition is only satisfied at I > 10 4 . 

In the limit of large sky coverage, one can Taylor ex- 
pand G2 around /, / =0, keeping terms up to f 2 , f ' 2 . 
The resulting 2-point joint PDF takes the Gaussian form 



Pi ex exp 



X 2 + X 



2XX r 



2(1 -r 2 ) 



where X = 



C-C 
o-c 



(38) 

and r is given by equation ([5]). If £ ~ t, then Gg and 
Cgi are dominated by contributions from halos with sim- 
ilar M, z and Xg. So we expect Cg and Cy to be highly 
correlated. Their correlation is quantified by r and as ex- 
pected, when £ — > £, r — ► 1. When i differs significantly 
from £, Cg and C e < are dominated by different clusters, 
so r — > 0, meaning G/ and Cf are nearly independent. 

3.6. Numerical evaluation of the PDF 

The integral form of the n-point SZ power spectrum 
PDF (equation |3"U|) must be computed numerically. This 
is straightforward for the 1- and 2-point pdfs, but more 
sophisticated Monte-Carlo methods must be used for ef- 
ficient evaluation when n is large. In what follows, we 
illustrate our results using the 1-point pdf ( equation I23p . 

The solid curves in Figs. [H [5] and [6] show equation (|23[) 
for surveys covering 1, 10, 200 and 4000 deg 2 . For com- 
parison, we also show the corresponding Gaussian form 
( equation [36]) . For 1 deg 2 sky coverage, representative 
of CBI, P(G) is highly non-Gaussian at £ < 4000. This 
is consistent with what would be expected from Fig. [3l 
which shows that a significant fraction of the relevant 
modes deviate from the Re Yl Dj ^ f 2 scaling, so the 
Gaussian condition (equation l37j) is not satisfied. Fur- 
ther understanding of this non-Gaussian behavior comes 
from Fig. [TJ the effective number of clusters contributing 
to the SZ effect is N P / sky < 10 for £ < 4000. So we ex- 
pect large Poisson fluctuations in the SZ power spectrum 
for £ < 4000. Although N p increases with I, making the 
non-Gaussianity weaker at higher £, it is still only ~ 30 
at £ — 10 4 , so we expect P(C) to be mildly non-Gaussian 
even at I — 10 4 . For 10 deg 2 sky coverage, roughly the 
area SZA 3 plans to cover, P(C) is strongly skewed at 

3 SZA, http://astro.uchicago.edu/sza/ 
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Fig. 4.— The PDF of C = Q for £ = 1085,2270,4745,9923, 
respectively in a 1 deg 2 survey. Solid lines show the pdf given 
by equation H23I I, and dashed lines show a Gaussian form with the 
same mean C and ctq for comparison. P(C) is highly non-Gaussian 
for I < 4000, mainly because the small survey volume means that 
discreteness fluctuations in the cluster counts are large (c.f. Fig. 
0. 



/ < 2000, with a non- negligible tail at high C. Only for 
multipoles > 5000, where N p f sky > 100 (Fig. Q}, does 
P{C) approach the Gaussian form (Fig. [5]). For larger 
sky coverage, P(C) becomes more Gaussian. For a 200 
deg 2 survey, roughly ACT 4 and APEX 5 plans to cover, 
P{C) approaches Gaussian at I > 2000. Note, however, 
analysis of the SZ power spectrum at I < 2000 must ac- 
count for the fact that the pdf is non-Gaussian. Again, 
these results are consistent with expectations based on 
Figs. [T]and H 

For 4000 deg 2 sky coverage, roughly the size of SPT 6 , 
P(C) is Gaussian at £ > 1000. For Planck, which will 
cover the full sky, P(C) can be assumed to be Gaussian 
distributed at all relevant scales. While this will signif- 
icantly simplify the analysis of the SPT and Planck SZ 
power spectra, we caution that the variance of this Gaus- 
sian differs significantly from that expected in a truly 
Gaussian random field (fig. [1]). The strong correlation 
between different £ modes makes the total S/N of the 
SZ power spectrum significantly smaller than that of a 
random Gaussian field. 

In our model, the SZ power spectrum is dominanted 
by contributions from single halos (the one-halo term), 
and we have attributed the non-Gaussianity in the SZ 
p(Cg) we find to Poisson fluctuations in the halo number 
distribution. In this model, different I modes of the SZ 
effect can be highly correlated, with correlation length 
A£ c comparable to £ (Fig. [2J. Hence, the PDF of the 



4 ACT http:/ /www. physics. princeton.edu/act/ 

5 APEX. http:/ /bolo. berkeley.edu/apexsz/indcx. html 

6 SPT, |http:/ /spt.uchicago.edu/ 1 



Fig. 5. — Similar to Fig. [4] but for 10 deg 2 sky coverage and £ = 
520, 1085, 2270, 4745, respectively. The non-Gaussianity is weaker 
compared to that in Fig. f3] because the larger number of clusters 
in the survey volume brings the Poisson fluctuations closer to the 
Gaussian limit. 



band power with width A£ is roughly the same as the 
PDF of a single mode, as long as A£ < A£ c , The strength 
of the non-Gaussianity in this band power PDF is then 
also determined by the effective number of halos in the 
survey volume, which is iVp/ s ky. 

However, there is another possible origin of the non- 
Gaussian PDF. Suppose that the distribution of SZ tem- 
perature fluctuations were exactly Gaussian. Then p(Cg) 
would follow a ^-distribution, so one might wonder if 
this is the origin of the non-Gaussianity seen in Figs. HJ 
[5]and[6]? The argument for band-powers, rather than sin- 
gle modes, is similar: If Yi is the i-th Fourier mode, and 
we define the normalized band power x = |i^| 2 /(nC), 
where C is the mean variance of the Yi, then x is the av- 
erage of n independent Fourier modes. If we approximate 
the variance in each mode as being the same as C, then 



P n (x) = 



■x n ~ x exp(-2;/2)2n 



(39) 



is a x 2 -distribution with In degrees of freedom (recall 
each Yi has independent real and imaginary parts). Note 
that P n (x) becomes increasingly Gaussian in shape as 
n ^> 1, but it is very non-Gaussian for small n. 

For a survey with fractional sky coverage / s ky, n — 
£A£f s ky, where setting the bin size Al equal to the 
correlation scale A£ c is a natural choice. Hence, if 
n — £ A£ c / s ky 3> N p j ' s ky, then the non-Gaussianity in 
the band power PDF will be mainly contributed by Pois- 
son fluctuations in halo number. Now, Fig. [5] shows that 
A£ c ~ £, and Fig. [T] shows that, at £ ~ 10 3 , 7V P ps 



WW 



so n 



£ 2 is indeed larger than N p w (£/Wy 



This shows that the non-Gaussian shape of p(Ce) in our 
model is dominated by the Poisson fluctuations in halo 
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Fig. 6.— Similar to Fig. \E\ but for 200 and 4000 deg 2 sky 
coverage and £ = 520, 1085, respectively. Comparison with the 
smaller surveys shown previously (1 and 10 deg 2 ) shows that the 
pdf becomes increasing Gaussian as the sky coverage increases. 



(cluster) counts; Fig. [7] shows this explicitly. 

The PDF of the SZ power spectrum is similar to that 
expected from gravitational lensing in several respects: 
(1) P{C) is skewed so that it peaks at C < C; (2) It 
has a long tail extending to high C. For BIMA/CBI, 
this increases the probability of having C > 2C at the 
relevant £ range. When combined with the large sample 
variance caused by strong correlations between different 
£ modes, this may help reduce the tension between the 
CBI/BIMA power excess relative to the naive expecta- 
tions based on primary CMB and the Gaussian model 
of P(C). We postpone a quantitative discussion of this 
issue to elsewhere. 

However, the angular scale dependence of the SZ power 
spectrum PDF is distinctively different from that due 
to gravitational lensing. The lensing Ce is dominated 
by the two-halo term at large scales, whereas the SZ 
power spectrum is always domina ted by the one-halo 
term (Ko matsu fc K itavamal [X999r ). Hence the pdf in- 
duced by lensing is expected to become Gaussian at large 
angular scales (small £), whereas the SZ PDF becomes 
Gauss ian at small angular scales (large £) (jZhang fc Pen! 
(2001) and see Figs. HI El and[(5]). Furthermore, the large 
scale SZ power is mainly contributed by massive clusters, 
whose number density is low (e.g. iKomatsu fc Seljakl 
((20021) and Fig. []} . This causes large Poisson fluctu- 
ations and strong deviation from Gaussianity at large 
scales. On smaller scales, where more and more small 
clusters contribute, the Poisson fluctuations tend to the 
Gaussian limit. 

A remaining question is the dependence of p(Ci) on 
cosmological parameters, especially as- A convenient 
measure of p{Cg) is the effective number of clusters iVp 
(equation [5]) . There are essentially two effects: (1) A 



Fig. 7. — The origin of the non-Guassianity in the SZ band 
power spectrum PDF. If the SZ signal were Gaussian distributed, 
then the PDF of the SZ band-power would follow a x 2 -distribution 
with 2£A£/ S ]jy degrees of freedom, and one might wonder if this was 
the origin of the non-Gaussian shape of the SZ PDF. To illustrate 
that this is not the origin of the non-Gaussianity in our model, the 
long-dashed lines show Chi-square distributions when (£, At) = 
(520,400) and (1085,800), respectively. They are clearly different 
from the solid curves which show the PDF predicted by our model; 
short dash lines show a Gaussian form with the same rms. (In fact, 
because close £ modes are strongly correlated, the true PDF of the 
SZ band power with bin size A£ < A£ c is roughly the same as the 
PDF of a single Fourier mode (Fig. l4l l5l and \Ei . so we have simply 
copied the corresponding curves here.) 



change in er 8 changes the concentrations of halos, and 
hence relative weights of clusters. Since Np is the 
number of clusters weighted by their SZ contributions, 
change in c affects iVp . (i) Roughly, a cluster contributes 
to the SZ power spectrum only at scales larger than 
r s = r vu-/c cx M 1 / 3 /c; the SZ signal from a cluster is 
smooth on scales smaller than r s . Here, c is the concen- 
tration parameter and r v j r is the virial radius of a clus- 
ter of mass M . Since c decreases with increasing mass, 
massive clusters only contribute to the SZ signal on rela- 
tively large scales; whereas lower mass clusters contribute 
down to smaller scales. Now, increasing as makes clus- 
ters of a fixed mass more concentrated — at fixed mass r s 
decreases. Thus, a massive cluster which does not con- 
tribute power to a given scale if as is low may be able 
to contribute at higher as- Thus, increasing as increases 
the contribution from rarer more massive clusters, (ii) 
Increasing as is exp ected to increase c by the same frac- 
tion for all masses dBullock et al.l 120011 : IWechsler et al.l 
12001 iZhao et al.l 12003). However, a constant fractional 
increase in c results in larger fractional increase in the 
SZ flux from a more massive cluster. Thus the weight- 
ing of more massive clusters in the SZ power spectrum 
is increased when increasing as- (2) Increasing as in- 
creases the number density of massive clusters, making 
A^p larger and so Poisson fluctuations in cluster abun- 
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Fig. 8. — Dependence of the effective number of clusters Np 
(equation [5J on erg. The short-dashed (erg = 0.8) and dot-dashed 
lines (erg = 0.7) have the same relation between concentration c 
and halo mass m as when erg = 0.9. Hence, the difference in Np is 
entirely due to the change in the mass function dn/dm. At large t, 
where less massive clusters dominate, the dependence on erg is neg- 
ligible. This is because the abundance of the most massive clusters 
is exponentially sensitive to <rg, whereas less massive clusters are 
less strongly dependent on erg. The dotted (erg = 0.8) and long- 
dashed (erg = 0.7) lines account for the additional fact that halo 
concentrations also dependend on erg . This dependence partly can- 
cels that from dn/dm. At small I, the effect on dn/dm dominates; 
at large t, the effect on c dominates. Hence, at I ~ 10 4 , Np when 
erg = 0.7 is larger than when erg = 0.9. 



dances smaller. 

Effect (1) means that increasing as increases the rela- 
tive contribution of the rarer more massive clusters (by 
increasing c) ; if cluster abundances were not also altered, 
this would decrease Np (since massive clusters are rarer). 
However, increasing as also increases the number of more 
massive clusters (effect 2), so the two effects act in ap- 
proximately opposite directions. Figure [8] shows the re- 
sult of an explicit calculation of how Np depends on as- 
the two effects do indeed approximately cancel. 

When expressed as a function of x = Ci/Ce, the 
PDF p{Cg) is approximately invariant to changes in as- 
This provides some a nalytic support for the procedure of 
iDawson et all (|2006l ) which we discussed in the Introduc- 
tion. However, our analysis indicates that this procedure 
is only approximately correct. In reality, effects (1) and 
(2) do not cancel exactly; we do see differences in N p and 
p(Ce) as as is varied. Robust evaluation of the net effect 
requires improvements to the models we use for cluster 
abundance, the c — M relation and the gas model, as 
well as the incorporation of several complexities which 
we discuss in the next section. These improvements will 
be addressed elsewhere. 

4. SUMMARY, DISCUSStON OF LfMITATfONS, AND 
EXTENSIONS 



We have derived an analytic estimate of the n-point 
PDF of the SZ power spectrum (equation [30]). Our 
derivation is based on the halo-model. The analytical 
integral form of the PDF allows fast calculation of the 
one- and two-point pdfs, P(Ci) and P(Ce,Ce>) (equa- 
tions [23] and l27f . More advanced integration routines 
are required to compute higher-order PDFs in reason- 
able timeframes. 

We find that the non-Gaussianity of the SZ power 
spectrum is a function of scale i, and, at fixed is a 
strong function of the survey area (Figs. 01 [5] and [6J. 
This strong non-Gaussianity may resolve the discrepancy 
between CBI/BIMA SZ measurements and the primary 
CMB measurements. For survey areas ~ 10 deg 2 , which 
will be achieved by ongoing surveys such as SZA, the 
non-Gaussianity is significant at I < 4000, but becomes 
negligible at I > 5000. For future SZ surveys such as 
SPT, which will cover 10% of the sky, the PDF of the 
SZ power spectrum can be approximated as Gaussian at 
£ > 1000. For Planck, the PDF of the SZ power spec- 
trum can be approximated as Gaussian at all relevant 
scales. 

While useful, our results should be treated with cau- 
tion for a number of reasons. (1) We have assumed clus- 
ter dark matter density profiles are universal, and can be 
described by the NFW profile. N-body simulations show 
that a non-ne gligible fraction of clusters are not well fit 
by this form rfjing fc Sutol 12000). (2) At fixed mass, ha- 
los which are w ell fit by the NFW form have a rang e 
of concetrations (jBullock et al.ll200ll ; Ijmg fc Sutdl2002l ). 
We have not accounted for this scatter. (3) We have 
assumed that clusters are spherical. In re ality, clusters 
are b etter described as tri-axial spheroids (jJing fc Sutol 
2002). This will cause the SZ s ignal of a cluster to be 
non-spherical ijWang fc Fan! |2004| ) . These complexities 
will increase the variance of the distribution of Cg. 

Of perhaps more concern is that fact that our assump- 
tion of a spherically symmetric profile leads to the un- 
realistic prediction that different £ with the same ampli- 
tude \£\ are completely correlated. On the other hand, 
the measured quantity is almost always a band-power 
averaged SZ power spectrum, where the average is over 
different I, so the assumption of spherical profiles should 
be accurate to lowest order. Nevertheless, it might be 
interesting to account for halo triaxiality. 

Incorporating (1), (2) and (3) into our model 
(equations l23l a nd l30|) is relatively straightforward. 
ICoorav fc Shethl ([2002) shows how the scatter in halo 
conc entrations can be incor p orated into the halo model, 
and ISmith. Watts fc Shethl (|2006f ) show how to incor- 
porate the effects of halo triaxiality into models of the 
dark matter power spectrum. Since this is particularly 
straightforward for the one-halo term, it might be inter- 
esting to combine their analyses with ours. But doing so 
is beyond the scope of our present work. 

While P(C) is important for assigning realistic errors 
to measurements of the SZ power spectrum, thus permit- 
ting accurate constraints on cosmological paramteres, it 
also contains important information about cosmological 
and gas physics. For example, the shape of the tail at 
C > C is sensitive to the number of clusters in the survey 
and to the gas fraction in clusters. With sufficiently large 
sky coverage, one can divide the survey area into many 
1 deg 2 patches, thus providing a direct measurement of 
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the SZ PDF on 1 deg 2 scales. Comparison with Fig. [4] 
then constrains cluster abundances and gas physics. A 
detailed analysis of how much additional information (i.e. 
more than is provided by the power spectrum itself) is 
contained in the full PDF is beyond the scope of this pa- 
per, although it certainly deserves further investigation. 

It would have been nice to test our predictions with 
measurements from numerical simulations. However, 
current SZ simulations lack sufficiently large independent 
simulation volumes to measure the PDF reliably beyond 
the lowest moments. That said, we stress that the vari- 
ance of the SZ power spectrum which our approach pre- 
dicts is in excellent agreement with simulations (c.f. $2]). 
This gives us confidence that our approach should at least 
provide a useful template for more accurate model of the 
SZ PDF. 

The analytical approach described in this paper can 
be applied to several other cases. For example, it can be 
applied to calculate the PDF of the lensing power spec- 
trum at high I straightforwardly, where the one halo term 



dominates. At low I where two halo term becomes im- 
portant or even dominant in the lensing power spectrum, 
this is less straightforward. Extending this approach to 
include the two halo term and so make predictions for 
the PDF of the weak lensing power spectrum is under 
investigation. 
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APPENDIX 
SMALL SKY COVERAGE LIMIT 
For surveys with very small sky coverage, | exp[J [e % ' c — X)dN\\ <C 1. One can then Taylor expand 

exp J (e^ c - l)dN ~l+f (e^ - l)dN + H^H^^l . 
Inserting into equation (|17p and Fourier transforming to real space gives P{C) in the Poisson limit: 

P{C) = (1 - N)5 D {C) + J 5 D (Cf sky - S)dN 

+ \ (/ Sd{Si + S2 ~ C ^y) dN ^ 

-2[ J 6 D {Cf sky - S)dN]N + N 2 S D {C) 

+ T (/ 6o( - Sl + S2 ~ C f^y) b i b 2dNidN 2 

-2[ / 6 D (Cf sky - S)bdN] [bdN+[[ bdN]' A S D (C) 



Here, N = J dN. Note that this expression does not rely on the assumption that Ps is Gaussian. One can Taylor 
expand equation (|18p in the limit / s k y — > 0, perform the Fourier Transform, and verify that equation (|18[) does indeed 
agree with the above result to order 0(a 2 ). This suggests that equation (fl~8f should work well for all sky coverage, 
even though the assumption of Gaussian Ps breaks in the small sky coverage limit. The Poisson limit requires 

Re(G)/ sky ~ Re(^ D s ) f sky « 1 . 



Fig. shows that only for / s k y <C 10 5 , or the sky coverage <C 0.4 deg 2 , will the Poisson limit be reached. In reality, 
the above Poisson limit expression of P(C) has very limited application, because all ongoing and upcoming SZ surveys 

sky 



have / sk y > 10 5 



